On surface-symmetric spacetimes with collisionless 
and charged matter 

Sophonie Blaise Tchapnda 

Abstract 

Some future global properties of cosmological solutions for the Einstein- 
Vlasov-Maxwell system with surface symmetry are presented. Global exis- 
tence is proved, the homogeneous spacetimes are future complete for causal 
trajectories, and the same is true for inhomogeneous plane-symmetric solu- 
tions with small initial data. In the latter case some decay properties are also 
obtained at late times. Similar but slightly weaker results hold for hyperbolic 
symmetry. 

1 Introduction and main results 

In general relativity the time evolution of self-gravitating collisionless particles 
can be modelled by the Einstein- Vlasov system. For surveys of results on that 
system see [T]. Cosmological spacetimes are those admitting a compact Cauchy 
hypersurface. In this case the particles are galaxies or even clusters of galaxies. 

Results on expanding cosmological models with collisionless matter, a posi- 
tive cosmological constant and surface symmetry have been obtained in |121 113] . In 
the present paper we want to examine what happens when the particles in the self- 
gravitating collisionless gas under consideration are charged. The Einstein- Vlasov 
system is then coupled to the Maxwell equations determining the electromagnetic 
field created by the charged particles. Considering this particular problem extends 
the knowledge on global dynamical properties of solutions of the Einstein equa- 
tions. Adding the Maxwell equations could also help to answer the question why 
in the cosmology literature people usually talk about the magnetic field more than 
the electric field. If the accelerated expansion, which here is due to the positive 
cosmological constant, could let the latter decay faster than the former then this 
could be an explanation. This might also have some connection with the so-called 
Landau damping effect [TS]. The results of the present paper do not suffice to 
address this issue, but we hope they can provide a basis for doing that in the 
future. 

As known results related to the Einstein- Vlasov-Maxwell system we can men- 
tion a small data global existence theorem in the spherically symmetric asymptot- 
ically flat setting obtained in [5] . In the absence of Vlasov matter the asymptotic 



behaviour of solutions of that system with T^-Gowdy symmetry was studied re- 
cently in [TT] . 

Let us now formulate our system. We suppose there are two species of charged 
particles, one of positive charge +1 and the other of negative charge —1. All the 
particles are supposed to have the same rest mass equal to 1, and to move forward 
in time so that the number densities and /" for positive and negative charge 
species respectively, are non-negative functions supported on the mass shell 

PM {ffa^pV = -1, P° > 0}, 

a submanifold of the tangent bundle TM of the space-time manifold M with metric 
g of signature — +++. We use coordinates (t, x"") with zero shift and corresponding 
canonical momenta ; Greek indices always run from to 3, and Latin ones from 
1 to 3. On the mass shell PM the variable becomes a function of the remaining 
variables {t,x°',p'') : 

P° = \/-9"°\/'^+9abP''p''- 
The Einstein- Vlasov-Maxwell system now reads 

dtf+ + ^9."/+ - ^(r^^PV + Fp V)9p./+ = (LI 



Gafi + ^Qap = ^T^{Taf3 + Tap) (1.3 

V^F^T, -f V;3F^a -I- V^F„/3 = (L4 
VaF"'^ = JP (L5 



Taf3^- / r +/ )PaPf3\g\ ' (L6 



= Fa^Fp ^ - i^F-'^F^s {1.7 



1/2 dp dp^dp'^ 
Po 



(1- 
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where Pet = ga^P^ , are the Christoffel symbols, \g\ denotes the determinant of 
the metric g, Gap the Einstein tensor, A the cosmological constant, F the electro- 
magnetic field created by the charged particles, the total particle current den- 
sity generated by the charged particles and Tap and Tap are the energy-momentum 
tensor for Vlasov and Maxwell matter respectively. 

A computation in normal coordinates shows that Vq, J" = 0. This equation 
is an expression of the conservation of charge. It can be shown as in [10] that T^p 
satisfies the dominant energy condition i.e. TapV°'W^ > for any two future- 
pointing timelike vectors and W"'. Let us show that the same is true for the 
Maxwell tensor Tap- Proving this is equivalent to show the weak-energy condition 
TapV^V^ > for all timelike vector V", together with the property that TapV^ is 
non-spacelike for any future-pointing timelike vector . The proof of the latter 
can be deduced from the following identities which hold since F is antisymmetric 

Tai^Tp ^'^{T'^^T^s)gap, TapT°''^ > Q 

Contracting the first of these identities twice with V"' implies the following, using 
the second identity and the fact that V"" is timelike : 

[V'^ra.){T''pV^) = \{T"T,s)gapV'^V^ < 

and setting Pi, — V^Tai^, this means that PaP" < 0, that is Pa is non-spacelike 
as desired. Now proving the weak-energy condition is equivalent to show that Too 
is non-negative since we can choose an orthonormal frame such that V" is the 
timelike vector of the frame. In such a frame goo — ~1 so that tqo = ^g°'''FoaFob + 
jF°-^Fab > as the sum of spatial lengths of a vector and a tensor respectively. 

In the present paper we adopt the definition of spacetimes with surface sym- 
metry, i.e., spherical, plane or hyperbolic symmetry given in [8j. We write the 
system in areal coordinates, i.e. coordinates are chosen such that R — t, where 
R is the area radius function on a surface of symmetry. The circumstances under 
which coordinates of this type exist are discussed in [2] for the Einstein- Vlasov 
system with vanishing A, and in [T3] for the case with A. The analysis there can 
be extended to the situation under consideration here since the Maxwell tensor 
Tap satisfies the dominant energy condition. In such coordinates the metric takes 
the form 

where 

r sin 6* if fc 1 

siufc 9 := < 1 if fc = 

[ sinh 6* if fc = -1 

Here t > 0, the functions A and fi are periodic in r with period 1. It can be 
shown as in [6] and [2j that due to the symmetry and /~ can be written as a 
function of 

t, r, w := e^p^ and L :== t^ip'^f + sin^. 9{p^f, with r, w e M ; L e [0, +oo[. 
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In these variables we have p° ~ e~''-\/l + + Lji"^. 

In surface symmetry the only non-zero components of F are Fq\ and F^-i- 
Indeed setting h :— g + ® — e\ ® e\, with cq = e^''^ and e\ — e^'*'^, the 
mapping X^^ ^ h^X°' is the orthogonal projection on the tangent space of the 
orbit, and since the vector Y„ := Faf3{eQ)°'h^ is invariant under the symmetry 
group, it vanishes. This implies that Fq2 = ^03 — 0. Similarly, replacing cq by ei 
in the expression of yields F12 = i^ia = 0. 

Now we can calculate the Maxwell equations in a coordinate frame. Equation 
(|1.4p then implies the following, where (9o, 01,82,83) = {8t, 8r, 8g, 8^) : 

80F23 = 81F23 = 82FQ1 = 83FQ1 = 0. (1.10) 

Using the fact that the mapping {p^,p^,p^) 1— > {w, L,p^) is one-to-one from IRxjO, oo[xI 
to RxjO, oo[xR and from Rx] — oo,0[x]R to Mx]0, cx3[xR, one can compute the 
J"'s and obtain = J3 = 0, 

/ {f + - .r){t,r,w,L)dLdw 

-CX3 Jo 



and 



^' = -Je-M / . : ^ (f+--nit,r,w,L)dLdw. 
Equation (|1.5p then implies 



d2{V\9\F^^) = 0, 83iV\g\F'^) = 0. 



(1.11) 



The non-zero components of the electric and magnetic parts of F are = e^F^^ 
and Bi = e~^-\/[gjF^^ respectively. Using these identities and recalling that 
— t^e^^^' siufc 9, equations ()1.10|) and (jl.lip lead to the following, where 
c is an arbitrary constant and is denoted by E: 

/OO 1*00 
I {f^ - r){t,r,w,L)dLdw, 
-00 Jo 

/"OC />OC 

dt{t'e^E{t,r))^-e^^n / -^=^=={f+ - r){t,r,w, L)dLdw, 

^1 + + L/i^ 



Bi(i,r) 



f-2pA(t,r) 



After calculating the Vlasov equations in the variables {t,r,w, L), the non- 
trivial components of the Einstein tensor, and the energy-momentum tensor and 
denoting by an upper dot or by prime the derivation with respect to i or r respec- 
tively, the complete Einstein- Vlasov-Maxwell system then reads as follows 

oA'-A, 



8tf+ + , --8rf+ - {\w + e>'-^^l'^Jl + w^TW- e^+^E)d^f+ = 

(1.12) 
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dtf- + ^====drr - (Aw; + e'*- Wl + + L/t^ + e>^+^E)d^f- = 

(1.13) 

e-'^i^{2t\ + l) + k- Kt^ = S,TTt^p (1.14) 
e-^''(2i/i-l)-A; + Ai2 = 87ri^p (1.15) 
^' = -A-Kte^+i^j (1.16) 

e-'^ (m" + - A')) - e-'" (^A + (A - A)(A + 1)^ + A = 47rg (1.17) 

drife^E)^fe^a (1.18) 
dt{t^e^E) = -t'^e% (1.19) 

where 

/OO /"CxD 
^Jl + ^JFTLj^{f+ +r){t,r,w,L)dLdw 

+ l(e2A^2 ^ct-4) = e-2''(Too +roo)(i,r), (1.20) 



/-OO ^oo 2 

P{t,r):=-^ / ^=====(/++/-)(t,r,«;,L)(iL(i«; 

- I(e2^^2 ^ct-4) = e-2A(rn +rn)(t,r), (1.21) 



j{t,r):=^j w{t + r){t,r,w,L)dLdw = -e^+^Toi{t,r), (1.22) 

/'OO /-OO 7- 

q{t,r):=-^ / ^=_==(/+ + /" )(t, r, i)dLd^i; 

+ (e^^E^ + ct-^) = ^{T22 + T22)(i, r), (1.23) 

/OO />00 
/ {f+-f-){t,r,w,L)dLdw, (1.24) 
-OO Jo 

(.OO (.OO 

b{t,r):=-^ / ^====(/+-/-)(t,r,«;,L)dLd«;. (1.25) 
We prescribe initial data at some time t = to > 0, 

o o 

f+{to,r,w,L) = f+{r,w,L), f-{to,r,w,L) = f-{r,w,L), 

X{to,r) = A(r), ii{to,r) = ti{r), E{to,r) = E{r) 

and want to study the existence and behaviour of the corresponding solution for 
t e [to, +oo). 
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To this end we maintain the notation in [71 [TU [T3] , and follow their work 
wherever possible. The first step consists on generalizing the local existence result 
in [71 theorem 3.1] to the case of charged particles under study: 

o o o 

Theorem 1.1 Let /± e C\M.'^ x [0,oo[) with /±(r + l,w,L) = f^{r,w,L) for 

o 

(r, w, L) £ X [0, oo[, > 0, and wq '■= Wq + Wq . Lq := L'^ + Lq with 

o 

:— sup{|w||(r, w, L) e supp/^} < oo 

o 

Lq := swp{L\{r,w, L) G supp/*} < oo 
Let °X, E e Ci(R), li e C^{R) with A(r) = A(r + 1), ^{r) = ^{r + 1), E{r) = 

o 

E{r + 1) for r e M, and 

l_i (r) = -47rtoe^+^ j(r) = e^+^ / / w{f+ + f-){r, w, L)dLdw, r e M, 

to J^oo Jo 

o Q o o fOO /'OC O O 

dr{tle^E) = tleK = 7re^ / {f^ - f~'){r, w, L)dLdw, r G R. 

J-oo Jo 

T/ien there exists a unique, right maximal, regular solution (f^ , f^ , X, fi, E) of 

o o O Q o 

jJ. jgp -fOgj) wit/i (/+, \,n,E){to) = {f^,f^,X,n,E) on a time interval [icTmaxl 

The regularity concept used in this statement is defined below. 

The next theorem provides a continuation criterion used to prove that the 
solution exists on the whole time interval [io,oo[. 

Theorem 1.2 Let [f^ , f^ ,\,^,E) he a right maximal regular solution obtained 
in Theorem If one has 

svcp{\w\\{r , w , L) e supp/^} < oo, sup{|i(j||(r, w, L) e supp/^} < oo, 

sup{/i(i,r)|r G K, i G [to,T,r^ax[} < oo, sup{(e^|£:|)(i, r)|r G R, i G [to,T„,ax[} < oo 

then Tmax = oo. 

In the following we claim that among the conditions given in the previous theorem 
there is one which implies the others. 

Proposition 1.3 The condition 

sup{/^(t,r)|r G R, i G [to,Trnax[} < oo 
is sufficient in order to conclude that Tmax = oo. 
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We can then state 

Theorem 1.4 Consider initial data as in Theorem \1.1\ and assume in the case of 
spherical symmetry that tg > 1/A. Then the solution of the surface- symmetric 
Einstein- Vlasov- Maxwell system with positive cosmological constant, written in 
areal coordinates, exists for all t G [to,oo[ where t denotes the area radius of the 
surfaces of symmetry of the induced spacetime. 

Once the existence of solutions is proven, one would like to study their asymp- 
totic behaviour at late times. In particular an important point is to know whether 
the spacetime obtained is future complete or not. This seems not to be easily 
achieved by a direct argument for the generic data case of the inhomogeneous 
Einstein- Vlasov-Maxwell system. Nevertheless it works out in the spatially ho- 
mogeneous case, as well as in the inhomogeneous plane-symmetric case under an 
additional assumption on the initial data. 

Let us examine spatially homogeneous solutions. These correspond to LRS 
(locally rotationally symmetric) models of Bianchi type I and type III and Kantowski- 
Sachs type for plane, hyperbolic and spherical symmetry respectively. We refer to 
|Dj for a detailed discussion on these models. We use the same notation as in [3], 
where the case of uncharged particles has been studied. Most of the results ob- 
tained in that case apply also if the particles are charged, the exception being those 
for which the proof involves matter terms. We will also use results from [14 in 
which the only requirement for the energy-momentum tensor is to satisfy the dom- 
inant and strong energy conditions, these are valid in the case under investigation. 
Let us formulate our system in Bianchi symmetry. 

The spacetime is considered as a manifold G x I, I being an open interval 
and G a simply connected three-dimensional Lie group. The metric has the form 

ds^ = — dr^ + gije^ (g) e-*, 

where {ci} is a left invariant frame and {e*} the dual coframe. 
The Einstein constraint equations are 

R - hjk'' + [hjg'^f = 16^(roo + too) + 2A 

= -SttToj. 

The evolution equations are 

dthj = -Rzj+(fci„.g'™)% - 2fcj;fc]- - 87r(r„- -I- r„) 

- 47r(roo + Too) + 47r(Ti™ -I- Ti^)g^"'g^, - Ag, 



(1.26) 
(1.27) 

(1.28) 
(1.29) 
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where 

Too + TOO = / (/+ + r)ir,v){l+9rsv'-v'y/^\g\'^'dv 



+ Fo^Fo -> + ^F^sF''^ (1.30) 

{f+ + r){r,v)vM^'dv (1.31) 

T,j+nj = J {f+ + f-){T,v)v,v,{l + grsv'-vT'^^lgl'^^dv 

+ F.^Fj T - ^F^sF-'\ (1.32) 

with V :— {v^ ,v'^ ,v^) and dv :— dv^dv'^dv^. 
The Vlasov equations are 

drf+ + [2fc>^- - (1 + grsv''v^)-"^f„,^v"^v- ~ {Fo + F, *4)]a„./+ = 

(1.33) 

5./- + [2A:>-'" - (1 + grsV^vn-'/^YmnV"'v- + iFo' + F, '^)Wr = 0, 



(1.34) 



where 7^„ = ^g'^ {-Cl^kdmi + C'L>5ni + Cl^ndki), C'jk are the structure constants 
of the Lie algebra of G. 

Note that in LRS Bianchi symmetry the only non-zero components of F 
are F^^ and F^^. The Maxwell equations (|1.4p allow us to obtain an explicit 
expression for the magnetic part F^^ of F. Thus the remaining unknown for the 
Maxwell equations is the electric part F"^ of F. For a Bianchi model the Maxwell 
equations ()1.5|) take the following form 

dA\9\'^'F'^]+Cl^F^0\g\'/^ = A9\'^' 

which yields for /3 = and l3 — i respectively 

Cl^F^° = J° (1.35) 

drF°' - itrk)F°' + CjijF^' = J\ (1.36) 

here 

J'^ f{f^-r)ir,v)\g\'/'dv 



r =/(/+- /-)(r, vyil + grsV^-vT'^'lgl'/'dv. 
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Setting kij = \{kimg^"^)gij + (Jij, crtj being the trace free part of the second 
fundamental form kij, and using the Hamihonian constraint (jl.26[) we obtain 

2 _R 1 

^ih^g'^f = — + ^(T,ja'' + 8np + A, 

then using the fact that R < (cf. [II]), it follows by the Proposition 1 in 0] that 
the matter energy density is bounded by 

p < Ce-^'-'^. (1.37) 

We can prove the following 

Theorem 1.5 Let f^{0,v) be a nonnegative function with compact support. 
Let (gij (0), /sij (0), /^(O, f ), /~(0, u), i^°'(0)) be an initial data set for the evolu- 
tion equations U.28\} . U.29]) . the Vlasov equations il.33]) . il.34\ l and the Maxwell 
equation U.36]) . which has Bianchi symmetry and satisfies the constraint equations 
lll.26\) . \L2'T^ , and \1.35\) . Then the corresponding solution of the Einstein- Vlasov- 
Maxwell system is a future complete spacetime for causal trajectories. 

In the inhomogeneous case the result in the latter theorem can be proved in plane 
symmetry for small initial data, a similar but slightly weaker result is true for 
hyperbolic symmetry as well. We have the following 

Theorem 1.6 Consider any solution of Einstein- Vlasov-Maxwell system with pos- 
itive cosmological constant in plane or hyperbolic symmetry written in areal coor- 
dinates, with initial data as in theorem \1.4\ 1^ the case of hyperbolic symmetry 
assume that A is sufficiently large. Let S be a positive constant and suppose the 
following inequalities hold: 

\toMh) -l\<d, |(e-V')(^o)| < -5, \{e^E){to)\ < 6 (1.38) 

|^^2g2p(to) _ 3 _ 3fce2*'(*«)| < 5, wito) <S, c < S, (1.39) 

where w{t) denotes the maximum of \w\ over the support of f~^(t) or f^{t). Then 
if 6 is sufficiently small, the following properties hold at late times: 

tX-1^ 0{t-^), e-^p' = 0{t-^), e^E = 0{t~^), (1.40) 

At^e'^'' - 3 - 3fce2^ = 0{t-^), w = 0{t-^). (1.41) 
Furthermore the spacetime is future complete for causal trajectories. 

The rest of the paper is organized as follows. In section 2 we present some pre- 
liminary results that we use to prove Theorem 1 1.1 1 in section 3. The proof for the 
other results is also given in section 3. 
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2 Preliminaries 

The regularity properties required for a solution are as in [7]. 
Definition 2.1 Let I c]0, oo[ be an interval 

(a) /± e C\I xM.^ X [0,00 [) IS regular, if f^{t,r + l,w, L) = f^{t,r,w,L) for 
{t,r,w,L) g / X X [0, oo[, > 0, and supp/^(i,r, ., .) is compact, uniformly 
in r and locally uniformly in t. 

(h) p (orp, j, q, a, b)e C^{IxR) is regular, if p{t,r+l) = p{t,r) for (i,r) e /xR 

(c) A e C^/ X M) IS regular, if X £ C^{I x M) and \{t,r + 1) = \{t,r) for 
{t,r) e I xR 

(d) p e C^(/ X M) is regular, if p' <E C'^{I x R) and p{t,r + 1) = p{t,r) for 
{t,r) e / X R. 

(e) E (or p) e C^{I x R) is regular, if E{t,r + 1) = E{t,r) for (t,r) G / x M. 

It is possible to solve each equation of the system p.l2p - p.l5p . (|1.19p . when the 
other unknowns are given. This is the content of the following 

Proposition 2.2 Let /+, f^, \, p, E be regular for {t,r) E L x R, I c]0,cxd[ an 
interval with to G /. Replace /+, f^, X, p, E respectively by /+, f^ , X, p, E in p, 

o o o ^ o 

b, p to define p, b and p. Suppose that f^ , f^ G C^(R^x [O, oo[), X, p, E £ C^(R) 
and are periodic of period 1 in r. Assume that 



t 

Then the system 



Stt A 



(2.1) 



dtf+ + , - {Xw + e^~^-p'^l + w^ + L/t^~ e^+f'E)d^f+ = 

^yl + W"^ + L/t'^ 

(2.2) 



dtf- + —====drf~ - {Xw + e'^-^^Vl + + L/t^ + e^+^E)d.^f- = 
•y/1 + + L/t'' 

(2.3) 

e~2A'(2tA + l) + fc-Ai2 = 87ri2p (2.4) 
e-2^(2t/i- 1) - fc + At^ = Stt^^p (2.5) 
dt{t^e^E) = -t^ef'b (2.6) 

o 

/las a unique, regular solution [f^,f^,X,p,E) on I x R with /+(to) — f^, 

_ °_ ° o ° 

f (^o) = / , A(to) ~ X, p{to) = p and E{to) = E. The solution is given by 

f^{t, r, w, L) = f^HR^, W^)ito, t, w, L),L) (2.7) 



10 



e 



t t ^ 3r 



A(t, r) = 4^ie2''(*''^)p(i, r) - ^^!^^| + -^te^f^<-''-^ (2.9) 

X{t,r) = X{r) + [ X{s,r)ds (2.10) 
J to 

£;(i, r) = t-2g-A(t,r) J^i2gl(r)|;(^) _ ^ e"'^' '%{s , r)ds^ (2.11) 

where {R^,W^) is the solution of the characteristic system 
d , . , ef"-^ 



^ {r, w) = (- ^ , , . - e^'^AVl + + L/t^ ± e^+f'E) (2.12) 

as ^yl + w^ + L/t^ 

satisfying {R^,W^){t,t, r,w,L) ~ {r,w). If I — [to,T[ with T G]to,oo] then there 
exists some T* e]to,T] such that condition 12. 1\) holds on [to,T*[xM. T* depends 
on p. 

Proof Integrating (|2.6p with respect to t over [to,t] gives p. lip . The rest of the 
proof is similar to those of propositions 2.2 and 2.3, 1) in [12 . □ 

In order to solve the system (|1.12p - (|1.19p . it will be enough to concentrate 
on the subsystem (|1.12l) - (|1.15p . p.l9p . as can be seen in 

Proposition 2.3 The subsystem U.12\ )- ri.l5\) . Iil.l9\) is equivalent to the full sys- 
tem il.l2]} - ![l.l9\} . provided the initial data satisfy 11.16\) and il.l8\) at t = to. 

Proof Under the assumption that the subsystem ()1.12p - ()1.15|) . p.l9p is satisfied 
as well as equations ()1.16p and p.lSp for t = tg, we should prove that (I1.16p - (ll.l8p 
hold for all t. Integrating (I1.19P over [to, t] with respect to t and differentiating the 
resulting equation with respect to r yields 

drit^e^E) = tle^a + ° s^g^^^' 4. ^'b)ds. (2.13) 

Using equations (|1.12p and (|1.13p and integration by parts with respect to t and 
w leads to an expression for s^e'^h'ds so that (|2.13p implies 

/•to /'OO /"OO 

dr{t^e^E)=t^e^a + Tr / / e^^+''E{f+ + f-){s,r,w, L)dLdwds. {2.U) 
Jt J -00 Jo 

Computing the conservation law V aJ" = in coordinates {t, r, w, L) and integrat- 
ing the resulting equation with respect to t yields 

/ / e^^+''E{f+ + f-){s,r,w,L)dLdwds ^0 

t J -oc Jo 
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and so ([^1^ becomes drife^E) = t^e^a that is PTT^ holds for all t. Using the 
latter and an argument similar to the one used in the proof for [7l Prop. 2.2] we 
can show that (|1.16p and p.l7p hold for all t as well. □ 

The latter proposition shows that equations (|1.16p and (|1.18p are invariant 
under evolution. So they will be considered as constraint equations on initial data 
(at t = to)- They can be solved : 

Proposition 2.4 The constraint equations ^ (r) = — 47rtoe '^^j(r) and 
drit^e'^E) = t^e^a are solvable. 

Proof To solve these equations we need to impose the following conditions, because 

o o o 

{e^E)[r) and e"^^'"' are periodic in r with period 1 : 

1 /"OO /"OO o o o 

A 



J-aoJO 



e (/+ - f'){r,w,L)dLdwdr = 0, 



and 



4_2 i-l j-oo pac o o 

I{f^):=— / / e^wf'^{r,w,L)dLdwdr = 

10 J-oo Jo 



to 

o 

Choosing A freely, the argument of the proof in ,23, Remark 2.4] applies. □ 

Remark Note that considering a model with more than one species of particles is 
important in order to prove the solvability of the second constraint equation above. 

1 oo oc ° ° 

Indeed if we had only one species of particles the integral e'^/(r, w, L)dLdwdr 

o 

would never vanish, except if / is identically zero. 

3 Proofs 

3.1 Proof of Theorem fTTl 

Instead of considering the subsystem (|1.12p - (|1.15p . (|1.19p . an idea used in [7], that 
we follow here, is to consider an auxiliary system consisting of 



dtf+ + , ^ .dr.f+ ~ {Xw + e^^-^n^l + w^+L/t^ - e^+'^E)d^f+ = 0, 

+ + L/t"^ 

(3.1) 



dtf- + -=jL^Z==drf- - {Xw + e^-^Vl + ^^'2 + L/t^ + e^+^E)d^f- = 0, 
^1 + + L/t^ 

(3.2) 

together with pTTI)) . pTTS)) . pTTg)l and 



fl = -iiTte^+^'j. (3.3) 
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Next by proving that fj.' — fl it is easy to show that if {f~^,f , X, ^, fl, E) is a 
regular solution of jSH]), pTM]) . pTTS)) . ((3:3)) . pTT9)) on some time interval 

/ C]0, oo[ with to S /, and with initial data satisfying (jl.lOp and (|1.18D for t ~ to, 
then {f+J-,\n,E) solves ^T^-^^. 

As in [7], the solution of the auxiliary system above is used to construct a 
sequence of iterative solutions. 

° o' ° o ° 

Let fL ^ , Ao(i,r) := A(r), ^loit,r) := ^(r), p.o{t,r) := fi{r), E(){t,r) := 

o 

E{r) for t G [fo,oo[, r G K ; Tq = oo. If A„„i, fln-i, En-i are already 

defined and regular on [to, r„_i[xR with T„_i > then let 

G„_i(i,r,w,L) := ( ==,-A„_iit; 



_ e^"-i-^"-i^„_i^l + w2 + ± e^"-i+^"-ii;„_i) (3.4) 

for t G [to, J^n^i [ and denote by (R^, Wjf){s, t, r, w, L) the solution of the charac- 
teristic system 

^(i?±,M/±) = Gti(s,i?,M/,i) 

with initial data 

{R^,W^){t,t,r,w,L) - (r,u;), {t,r,w,L) G [to,r„_i[xM2 x [0,c5o[. 
Define 

/±(t, r, w;, L) := /± ((i?±, T4^±)(to, t,r,w,L),L) , 
that is, is the solution of 

a/1 + 



Vl + w2 + L/t2 ^ e^-i+''"-ii;„_i)5^/± = (3.5) 



with /^(to) = and define p„, _p„, j„, 5„, a„, 6„ by the integrals (|1.20p -(|1.25p 
with E, A replaced by f^, -En-i, A„_i respectively. Define 

{ , , , ,to(e-2MW + Stt /■* 2 , N , 
T„:=sup-jt Gjto,7rj_i[| ^ k — J s pn[s,r)ds 



+ ^(i'-io) >0,rGM,tG [to,t']| 



Using Proposition [221 let 

to(e-2^('') +A:) , Stt /■* 2 , A 



e 



-2,„(,.) -0^- -k-^ s'p^{s,r)ds + ^(t' - 4), (3.6) 



t t ./to ' ' 3t' 
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1 -I- l-p'^t^ri(t,r) A 

A„(i,r) := 47rte2^"(*''-)p„(t,r) - + '-^te'>^-('^'-\ (3.7) 

Xn{t,r):=°X{r)+ f Xn{s,r)ds, (3.8) 

/2„(i,r) := -47rte^"+^" i„(t,r), (3.9) 

E„{t, r) t-2e-A.(t.'-) J^^SglM^^^^ _ ^* s2e'^"("''''6„(s, r)ds^ . (3.10) 

Subtracting equation p.Sp corresponding to by the one corresponding to 
integrating the resulting equation with respect to w and L and integrating by parts 
with respect to w yield 

drit^e^-bn) - (A„ - A„_i)e^"-i-^'-i+^"i2^„ - e^"-^-^"-^'-i+''"at(i2e^"a„), 

so that multiplying (|3.10|) by t^e^" and differentiating the resulting equation with 
respect to r lead to 



+ / s2(/i„-/i„_i)a„e^"-^'-i+^"-i ds. (3.11) 



J to 

We split the proof of theorem [Ll] into several lemmas. From now on || • || denotes 
the L°°-norm on the function space in question, the numerical constant C may 
change from line to line and does not depend on n or t or the initial data. Firstly 
we prove: 

Lemma 3.1 The sequences A„, A„, p„, Pn, En, j„, a„, hn, /i„e^"^^" are 
uniformly bounded in n, in the L°° -norm by a continuous function on [io,oo[. 

Proof Define P„(t) {P+ + P-){t) with 

P^{t) sup{|w||(r, w,L) G supp/,t(i)}, t G [io,T„[, and 

Qn{t) ■■= sup{se2^"("'''V e M,io < s < t}, 

5„(i) sup{|i;„|e^"(^^'')|r G M, to < s < t}, 
we have the following estimate on supp/,^(t) 



^\ + w^+L/t^< ^1 + (P±(i))2 + L±/t2 < c(l + Lo±)(l + P±(t)), 
so that 



II W IMI W II < + 4 )'(!+ II ll)(l + Pn{t)f{^ + Sn-l{t)f (3.12) 



II j^(t) II < ^(i + i±)^(i+|| /± ||)(l + P±(^))^ (3.13) 
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and then using (|3.7p and p.9p 

I e'^"-^>n(i,r) I + I K{t,r) \< C(l + A)(l + L±f (1+ || /± ||)(l + P±Wf 

(l + Q„W)(l + ^„-i(t))' (3.14) 

This inequahty is used to obtain an estimate on supp f^_^i{t) for iVF^^^j which 
impUes 

t 

to 



p.:t+i{t)<^o+c* I (i + p±(.))2(i + p±+i(5))(i + g„(.))3/2(i + ^„_i(,s)) 

) 



with = C(l + A)(l + L±)2(l+ II /± II). Setting P±(t) sup{P±(i)|m < n} 
and 5„(i) := sup{5m(t)|m <n}, it follows that 



P„^+i W < + C* (1 + P,f+i(s))3(l + Q„(s))3/2(i + 5„(s))3ds, 
whence 

Pn+iit) <wo + C* f {1 + P„+i(s))3(l + Q„(s))3/'(1 + 5„(s))3ds. (3.15) 
Taking the derivative of p.6p with respect to t leads to 

i(2A„e2'^") = 87r(te2'^") V + 4(te'^")' + yie^^" - A(ie2'^")2, 

integrating this over [to , t] and using integration by parts for the left hand side 
yields 

Qn{t) < ll^oe'^ll +C* (1 + P„(s))(l + Qn{s)f{l + Sn-i{s)fds. (3.16) 

Jto 

Next (IXTU)) implies that 

Sn{t) < \\e^E\\ +C* f {1 + P„(s))2(l + Qnis))'^^ds. (3.17) 

Jto 

Adding ([51^ . (PH)) and ([XTT]) implies that 

(Pn+i + 0„ + 5„)(i) < u;o + ||toe2^|| + ||e^l;||+C* / [l + {Pn+i + Qn + Sn)is)]'^^^ds. 

Jto 

Thus P„, Qm Sn are estimated by zi{t), the right maximal solution of the equation 
z,{t) ^wo + Wtoe^h + \\e^E\\ + C* f\l + z,{s)Y^'^ds, 
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to 



which exists on [to,^r'^^^[, we have T„ > T^^\ Therefore there exists a continuous 
function Ci(t) which depends only on zi as an increasing function such that 
K, K, Pn, Pn, En, jn, 6„, jlne>'"^-^" are bounded in the L°°-norm by Ci(t). 
□ 

Next we prove : 

Lemma 3.2 The sequences X[^, p'^, p'^, E'^, j'^, a'^, b'^, fi'^ are uniformly 
bounded in n, in the L°°-norm by a continuous function on [to,oo[. 

Proof Let us start from equation (|3.1ip . Given Lemma [5TT] most of the quantities 
on the right hand side of (|3.1ip can be estimated. The exception is ^i'^ which is 
obtained after differentiating p.6p with respect to r : 

fJ-'nit, r) = — ho/i (r)e-2^ - 47r s^p^s, r)ds 
again most of the terms are unproblematic with one exception that is p'^ 



oo />oo 



P'n{i, 7J / / y-, , "1 , r =7T^9r{f+ + /„ r, w, L)dLdw 

t J -co Jo y/l+W^ + L/t^ 

-e^"-'En-idr{e^"-'En-i). 

Defining 

Dnit) := sup{|| drf+is) II + II drf-is) \\ \to < s < t}, 

and 

A„(i) := sup{|| drie^"En){s) \\ \to < s < t}, 

we deduce the following estimates 

II iXW II <Ci(i)(A.W + A„_i(0) 
II Pnit) II <Ci(t)(ci+A.W + A„_i(t)), 

where ci :=|| e^^^/i || + || A || +1+ | A |, and so equation (j3.11|l implies 

\\dr{e^-En){t) \\<Ci{t)+ f Ci{s) {ci + Dnis) + An-iis)) ds 

JtQ 

thus ^ 

A„(i) < Ci{t) + f Ci{s) (ci +I?„(s) + A„„i(s))ds. (3.18) 

Jto 

Now differentiating (|3.7p and (|3.8p yields 

X'„it, r) = e2p„ (^STTtfi'^it, r)p^{t, r) + Antp'^it, r) - ^p'^it, r) + Mp'^ 

A'„(i,r)=A(r)+ / A;(s,r)ds, 

Jta 
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and using the expression for we obtain 

II tp'nit) II <Ci(t)(D„W+A„-l(t)) 

II AUt) II, II A;(t) II < Ci(t)(ci + D^(t) + A„_i(t)). 
On the other hand it follows from (|3.9p that 

and 

I {e^'--^-lL^)\t,r) |<Ci(t)(ci+A.W + A„_i(i)). 
We can now estimate the derivatives of with respect to r and w 

drG^{t,r,w,L) = iinn - A„)'e^"-^" 



^1 + W2 +L/<2 



dwG^{t,r,w,L) = (e 



l^n—X-n 



1 + 



(1 +1(;2 +i/i2)3/2' 



^1 + w2 +i/t2 



An), 



and thus 



I drGtit, r,w,L)\< Ci{t){ci + D,,{t) + A„_i(t) + A„(t)), 
\d^Gt{t,r,w,L) I <Gi(t), 

for t € [io,r'^'[, r e K, L e [0,Lo] and | w |< Differentiating the character- 

istic system with respect to r, we obtain 

W^+,){s, t, r, w, L) =drGt{s, Rt+„Wt+„ L).drRt+i{s, t, r, w, L) 
it follows that 



I —dr{Rt+i,Wt+,){s,t,r,w,L) I <Gi(s)(ci+i?„(s) + A„_i(s) + A„(s)) 

\driRt+i,W^+i){s,t,r,w,L) I, 
therefore by Gronwall's inequality we obtain, 

for (r, w, L) £ supp/+^i(t) U supp/+(t) U supp/,7+i(t) U supp/-(t) 



I dr {Rt+, , ){to,t,r,w,L)\< exp 



Gi(s)(ci + L'„(s) + A„_i(s) + A„(s))ds 
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The definition of implies that 

II ||<|| II sup{| driRt,W^){to,t,r,w,L) I \{r,w,L) e supp/+(t) U supp/" (i)}- 

Combining this with the previous inequahty and using the definition of Z?„ we 
obtain the following : 



Dn+iit) < (II II + II d(^r.w)f ll)exp 



Ci(s)(ci + L»„(s) + A„_i(s) + A„(s))ds 



(3.19) 



Let Dn{t) :— sup{Dm{t)\m < n} and A„(i) := sup{Am(i)|TO < n}. Then {Dn)r. 
and (A„)„ are increasing, therefore adding (|3.18|) and (|3.19|) implies 



^„+l(t) + A„(0 < + (1+ II d^r..^)f+ II + II II) 



exp 



Ci(s)(ci + i)„+i(s) + A„(s))ds 



Let Z2 be the right maximal solution of 



Z2{t) - Ci{t) + (1+ II II + II II) exp 



Ci(s)(ci + Z2(s))ds 



which exists on an interval [to,T'-^^C [to,T'-^H- Then it follows that 

-D„+i(t) + A„(t) < Z2(i), te[io,r(2)[, neN. 

Therefore there exists a continuous function C2 (i) which depends only on Z2 as an 
increasing function such that all the quantities estimated against I?„ and A„ are 
bounded in the i°°-norm by C2{t). □ 

The following lemma deals with convergence of iterates. 

Lemma 3.3 The sequences f+ , /,7, A„, /i„, E^, A„, /i„, /i„, p„, j„, a„, 6„ 
converge uniformly on every compact subset [tQ,T'^'^'^] C [io,2^^^''[ c"^ which the 
previous estimates hold. 

Proof Define for t e [io, T^^)] 

a„(i) := sup {II (/++! - /+)(s) II + II (/-+! - /-)(s) II + II (A„+i - A„)(5) || 



+ II {^ln+l - i^n){s) II + II (e-^"i;„ 



-ii;„_i)(s) II; to < s < t} 



and let C denote a constant which may depend on the functions zi and Z2 intro- 
duced previously. Then using Lemma |3. 11 we have 

II Pn+l{t) - pn{t) II, \\ pn+l{t) - Pn{t) ||, \\ Pn+l{t) - Pn{t) ||, 
II jn+l{t) - 3n{t) II, II a„+i(0 - an{t) \\, \\ fe„+i(0 - &„(t) ||< Ca„(i). 



18 



Using mean value theorem to estimate differences e-^"+i+^"+i — e'*'""'"'^" and e'^"+'^ 
gMn deduce from (13.71) and (13.91) that 



II A„+i(<) - Kit) II, II M„+i(i) - /i„(t) II < Ca„(i), 
p.Sp thus impUes 

II (A„+i - A„)(t) \\<C f an{s)ds. (3.20) 
By mean value theorem p.6p implies 

II (Ai„+i - /i„)(i) \\<C f an{s)ds. (3.21) 

"'to 

Reasoning as in step 3 for the proof of theorem 3.1 in [6] we can prove that 



\iR,W)++,-{R,W)+\ito,t,r,w,L), \{R,W)-^,-iR,W)-\ito,t,r,w, L) < C / a 

Jto 

which implies, using the fact that /jf was defined in terms of the characteristics, 
and by mean value theorem 



II if:+i - fzm II, II (/,:+! - fnm \\<c a^.,is)ds. (3.22) 

Jto 

From (|3.10p we deduce that 

II {e^-E,, - e^"-i£;„_i)(t) ||< C f a„_i(s)ds. (3.23) 

-'to 

Adding ((3:201) - ((3:23)) gives 

a„(t) <C (a„(s) + a„_i(s))ds, 

Jto 

then by Gromwall's inequality 

an{t) <C an-i{s)ds, 

Jto 

and by induction 

(-171 + 1 

an{t) < ^, for n e N, te [to,T^% 

We then conclude the uniform convergence of the iterates. In L°°-norm, A„ A, 

Hn ^ H, fln^ /i, /+, fn ^ E.D 

In order to prove that the latter limits are regular in the sense of Definition 
12.11 we need to show the uniform convergence of the derivatives of the iterates 
above. 
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Lemma 3.4 The sequences drf^, d^fn, drfn , ^^fn , ^'n^ Mn; An; dr{e^"^En) are 
uniformly convergent on [io,T^^^]. 

Proof Fix t e [to,T'^% \w\ <U, L < min(L+, Lg ), to < s < t. For d e {dr,d^} 
and s H-> {R^{s),W^{s)) the solution of the characteristic system associated to 
equation (|3.5[) in f^, define 

^tis) e(^"-^")(^'«^)5i?±(s,i,r,iz;,L), (3.24) 

vtis) ■■= dW^{s,t,r,w,L) 

+ (VTT^^^^TLT^e^"-^" A„)|(,_(^±_^±)(^_,_^_„_^))ai?±(s, t, r, w, L). 

(3.25) 

Reasoning as in step 4 for the proof of theorem 3.1 in [3] it can be proved with 
minor changes that for all e > 0, there is a non negative integer N such that for 
n> N, 

{\it+i~£.t\ + \vt+i-vt m<ce 

+ C f\\ Ct+i - I + I ^7,^+1 - vt \)ir)dT. 

Jto 

(3.26) 

By Gronwall's lemma it follows that the sequences and ry^ converge uniformly. 
The transformation from {dR^ , dW^) to (^^, r/^) being invertible with convergent 
coefficients, we deduce the convergence of dr,w{Rt ^^t) ^^"^ thus the convergence 
of drfn dwfn^ usiug the fact that /jf was defined in terms of characteristics. 
Let us prove the convergence of A^, /i'j, dr{e^" En). Define 

7„(t) sup{| e++i - e+ 1 1 - e,; i 1 ^++1 - I (s) 

+ I Vn+l - Vn I II (m'„+1 - KJ(s) II + II (K+l - II 

+ II dr{e^"+'En+i) - drie^"En) \\ (s),to <s<t}. (3.27) 

The sequences fin, An? Am A„, e'^"En, Pn, jn, cin, bn converge uniformly, we then 
take the above integer N large enough so that we have for n > N, 

II {e^"En - e^"-'En-l) (s) ||, || (fin+l - Pn)is) ||, || (An " An-l)(s) II, 
II (An+1 - An)(s) II, II {pn+1 - Pn){s) ||, || {jn+1 - jn){s) ||, 

II (An+1 - An)('S) lUI (fln+l " an)(s) ||, || (^n+l - &«)(«) ||< £• (3.28) 

Taking d = dr,it follows from ((3:241) - ((3:251) that 

5i?±(s):=e(^"-^")(^-«J)e(s), (3.29) 
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dWtis) = Vnis) - {Vl + ^' + L/s^X,,)^t{s)- (3.30) 

Using these equations and since was defined in terms of characteristics we 
obtain 

II (a./±+i - II <|| drj^ II {\drRt+^ - drRtl + \drWt+, " drW^\){s) 

<\\ drj^ II (|e^'-^-""+^Ct+i 

+ - vt\ + |A„+ie;^+i - Ut\){s) (3.31) 
this implies, using (I3.28P and the fact that A„, A„, ^„ are bounded: 

II idrf^+, - drft){s) \\<Ce + C(|C±+i - C,tl + - '7,tl)(s). (3.32) 
Using p.27p . (|3.28p . p.32p and the expressions of p„, p„, j„, we deduce that 

II ip'n+l - Pn)is) \\,\\ (K+l -Pn)(s) ||,|| ( Jri+l " ) («) II < C'e + C'(7„ + 7n- 1 ) (s) • 

(3.33) 

Now taking the derivative of p.6p with respect to r, using p.28p . (|3.33p and the 
fact that /i„, /^^ are bounded we deduce an estimate for /x^: 

U^^Ul~PnM\\<Ce + C [\ln+ln-l){T)dT. (3.34) 

Jto 

For AJ^, we first take the derivative of (|3.7p with respect to r and obtain 

a; = {8T:tfi'„pn + 47rtp;)e2'^" + {At - k/t)fi'„e^^'" , (3.35) 

this shows that A^ is bounded. Subtracting (|3.35p written for n + 1 and n, we 
obtain the following, using the fact that p„, /i„, fi'^ are bounded and p.27p . (|3.28p 
and 

II (K+i^KM ||<C6 + C(7„+7n-i)(s), (3.36) 
and integrating this over [tg, t] it follows from p.Sp that 

II (A'„+i-A;)(s) \\<Ce + C A7„ +7„_i )(r)dr. (3.37) 

Jto 

For /i^, we take the derivative of (|3.9p with respect to r, subtract the expressions 
written for n + 1 and n, and use (|3.27p . (|3.28p and (|3.32p to obtain 

II (AUi-/i;j(s) ||<^^e + C(7„ + 7n-i)(s). (3.38) 

Now for an estimate for dr{e^" En), we subtract the expressions written for n + 1 
and n from (|3.1ip . use (|3.27p . p.28p . p.32p and the fact that A„, p'^, an, bn, 
jln are bounded to obtain 

II (9,(e^"+ii?„+i)-a,(e^"i?„))(s) \\<Ce + C f {in + ln-i){T)dT. (3.39) 
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Now we can add p.26p , (|3.34p , (|3.37p and p.39p to obtain, after taking the supre- 
mum over s G [to j ^ 

ln{t) <Ce + C f (7„ + 7n-i)(s)ds, (3.40) 

Jto 

and setting 7„(t) = sup{7„i,TO < n}, we deduce by Gronwall's lemma that 

ln{t)<Ce, n>N, te[to,T^% 

Thus the sequence 7„ converges uniformly to 0. By (|3.34p . (|3.37p . (|3.38p and (|3.39p . 
the sequences /i^, /x^, and dr{e^"En) then converge uniformly on [tojT'^^^]. □ 
The regularity of f'^, f~, A, fi, jl and E follows. Using the convergence of the 
derivatives, it can be proven as in [B] that (/+, A, /i, /i, E) is a regular solution 
of (lO). (IX^ . (fTTTD . (frra . dOl . (irrgi) and then (/+,/-, a, ^,i;) solves (HHD- 
l|1.19p . To end the proof of Theorem 1 we prove the uniqueness of the solution. 
Let {f^,f,^,Xi,fii,Ei), I = 1,2 be two regular solutions of the Cauchy problem 

o o o ^ o 

for the same initial data (/^, A, /i, _B) at t = tQ. Setting 

ait) := sup {II (/+ - /+)(.) II + II (/r - f^is) II + II (Ai - A2)(.) || 
+ II il^i - ^i2)i■s) II + II [e^'E, - e^-E.2){s) ||; to < s < <} , 

and proceeding similarly as to prove the convergence of iterates leads to 

a{t) <C f a{s)ds, 

Jta 

which implies that a{t) — for t G [to, This proves uniqueness and completes 
the proof of Theorem 11.11 

3.2 Proof of Theorem [lH 

Let (/"•", /~, A, /i, i?) be a right maximal solution of the full system (I1.12p - (|1.19p 
with existence interval [to,T,„Qa:[- We assume that Tmax < oo. By assumption 

Q, := sup{te2^(*''^)|r e M, to < t < T^ax} < oo, 
S, := sup{|^|e^(*''^)|r e M, to < t < T„,ax} < oo, 

and < oo where P* := P+ + with 

:= sup{|w|, {r,w,L) G supp/='=(t), t € [to,T„a^[}. 

We take ti (z\to,Tmax[, and we will show that the system has a solution with 
initial data (/"'"(ti), /~(ti), A(ti), /z(ti), i?(ti)) prescribed at t = ti which exists 
on an interval [ti, ti + S] with 6 > independent of ti. By moving ti close enough 
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to Tmax this would extend the initial solution beyond T^ax, a contradiction to 
the initial solution being right maximal. We have proved previously that such a 
solution exists at least on the right maximal existence interval of the solutions zi, 

Z2 of 

Zlit) = W{h) + ||tie2^(*i)|| + \\e^^'^^E{h)\\ + C* f\l + z,{s)Y'''^ds, 



Z2{t) = C,{t) + (1+ II II + II |l)exp 



Ci(s)(ci + Z2{s))ds 



where W{ti) := W+{ti) + W-{ti), 

W^ih) sup{|?«|, (r, w, L) G supp/±(ii)}, 

C* = C{1 + A)(l + Lo)'(l + ||/+(ii)|| + lir (ii)ll), 
ci :=l + A+||A'(ii)|| + ||e-2Mti)^'(t^)||^ 

and Ci is an increasing function of zi. Now W{ti) < P^, Utie^^'^*^^ || < Q^, 

o 

||e^(*i)i;(ti)|| < S*,, ||/^(ti)|| = ||/=^||, Lo is unchanged. Thus we have uniform 
bounds W{ti) + ||iie2^(*i)|| + ||e^(*i)£;(ti)|| < Mi, C* < Mj. On the other hand 
we can use the expressions for A', /i', A', some estimates proved in lemmas 3.1, 
3.2, and 3.3 to obtain uniform bounds ci < M^, Ci{t) + (1+ || d(r.w)f^{ti) II + I 
d(r.w)f~{ti) II) ^ M4. Let yi and y2 be the right maximal solution of 



yiit) = Ml + M2 / (1 + 2/i(s))i5/2^s, 



y2it) = M4exp 



CU.s){M3 + y2is))ds 



respectively, where C* depends on j/i in the same way as Ci depends on zi. Then 
2/1 and y2 exist on an interval [ti,ti+S] with 6 > independent oiti. If we choose 
ti such that Tmax < ti + S then zi < yi, Z2 < y2, in particular zi and Z2 exist on 
[ti,ti + 5]. This completes the proof of theorem 11.21 



3.3 Proof of Proposition 1.3 

We start by showing how to obtain the bound on w. Since we are in the non- vacuum 
case one has > and > 0. For t > to define 

Pf{t) := max{0, max{w|(r, w, L) e supp/^(t)}}, 
P^{t) := min{0,min{w;|(r, w, L) e supp/±(t)}}. 
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Let {r{s),w{s), L) be a characteristic curve in the support of Assume that 
P+{t) > for some t G [to,Tmax[, and let w{t) > 0. Let ti e [to,t[ be defined 
minimal such that w{s) > for s € [ti,t[. We have 

w{s) = ~\w - e/'-^n'^/l + w2 + L/s2 ;p e^+^'i; 
47r2 



s 



OO />CXD 



1 



{k - As2)e2M 



-w 



(3.41) 



2s 2s 

As long as w{s) > we drop the last two terms of the right hand side in (|3.4ip 
since they are negative, and then obtain 

3±/ 



w{s) < — 
s 



C w{l + L) 



w 



(/+ + Ddldw + - + e^+^'\E\ 
s 



< 



C 



w{s) 



w{s) 



\E\. 



(3.42) 



Now integrating ()1.19|) with respect to t, using the boundcdncss of /i and the fact 
that s2|6(s)| < s'^d{s) < C'Pf{s) (a being similar to a with /+ — replaced by 
/+ + /^), we obtain 

(e'+^mis) < Ctjie^mih) + C r Ptir)dr. 



It then follows that 



C 



w{sr < -(Pf (s))^ + CPf (s) + CPf (s) Pf (r)dr, 



ds s 
which implies after integration 



'{pt{s)r+pt{.^)] ds+c 



P^s)P^T)dTds. 



tl Jti 



If tl = to then w{ti) < wq, otherwise w(ti) — 0. In any case it follows that 



w^{t) <wl+C / [s-i(P±(s))2 + P±(s)] ds + C 



Pf{s)Pf{T)dTds. 



tl Jtl 



(3.43) 



The double integral in the right hand side of (|3.43p needs to be worked out. 

f r Pt{s)Pt{T)dTds<\ f [\p^{s) fdTds + \ f l\pt{T)YdTds 
Jt^ Jt^ ^ Jtl Jtl ^ Jtl Jtl 

ft 



(s-ti)(P±(s))2ds+i^ {Pt{T)fdsdT 



< Ct / (Pf (s))^ds. 
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Therefore (|3.43p implies 

{Pt{t)f <{wl + Ct) + C / (s-i +t + l)(P±(s))2ds, forte [to, T„ 



By Gronwall's inequality it follows that is bounded on [io, T„iax[- 
Estimating w^s) from below in the case it;(s) < along the same lines shows that 
is bounded as well. 
The bounds on w and fji imply that e^\E\ is bounded as well, using (|1.19p . 

3.4 Proof of Theorem 11.41 

We prove that /i is bounded on [tQ,Tmax[- 
A lengthy computation leads to 



< 



2p + q-^[l + {k^Ke)e"^] 



dr 



l + {k- At2)e2^]dr, 
(3.44) 



since q is nonnegative. Using the fact that p + p>0, p>p and k — At^ < it 
follows that 



dt Jo 

and by Gronwall's inequality 



e''+^p(t, r)dr < 



1 



oM+A 



pdr, 



[ e^+^p(t,r)dr < Cr^ t e [to,T^axl 
Jo 



(3.45) 



On the other hand using the equation p' = —Ante'^'^^j, (|3.45p and the fact that 
IjI ^ P we find 



\p{t,r)- / p(t,a)da\ <C, te [tQ,Tma 
Jo 

Next using (|1.15p . p — p < and k — At^ < 0, we have 
d 



e [0,1]. 



(3.46) 



dt 



gM-A _ gM-A 



1 -I- i-p^M 



t 



< e 



,Ai-A 



1 k-At^ 



t 



,2m 



t 



- t 
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so that 

^i^-Xmr) <Ct, te [to, Tr^aA, T G [0, 1]. (3.47) 

Now using p — p < and k — At'^ < 0, it follows that 



1 /■! pt pi 

li{t,r)dr — / p{r)dr + I / (i{s,r)drds 
Jo Jta Jo 



<C+ [ — [ [e^''(8Trs^p + k~As^) + l]drds 
J to 2s Jo 

<C+ [ ^1 (Sns'^e^'-^e^'+^p+Udrds 
J to 2s Jo 



< C + i In(Vio) + C f f e^+^pdrds 
2 Jto Jo 

<C+^ln{t/to) + Ct^, 
and using (|3.46p we obtain 

p{t,r) <C(l+<2 + lnf), te [ta,Trnax[re [0,1]. 

p is then bounded on [to, Tmax[ and by proposition 1.3 the proof of theorem ll.4l is 
complete. 

3.5 Proof of Theorem [HI 

The equation of motion for charged particles is given by the following differential 
system for a path v ^ (r, , v^){v): 

av av av ■' 

For a particle with rest mass m moving forward in time, v'^ = (to^ + gijV^v^Y^"^ . 
Then the relation between coordinate time r and proper time v is 

^^[m'+g.^^y^y^f/2_ (3.43) 

In order to prove the completeness of trajectories it is useful to control gijV^v^ as 
a function of r. As in 4 , we can define, from the Vlasov equations (|1.33p - (|1.34p . 
the characteristic curve V^'(t) satisfying 
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In this section we use the notation 7 :— and C is an arbitrary positive 

constant which may change from line to line. 
Using the latter equation we obtain 

< (-27 + Ce-'^)g'^V,V, T ^g''V,Foj. (3.49) 

We have used [H (3.22)]. The second term in the right hand side of p.49p can 
be estimated using the Cauchy-Schwarz inequality and the elementary inequality 
xy < + y^/2e : 

g^'V^Fo, < {g'^V^Vjf'\g'^F^,Fo,f'^ 

< '^(g^'V^V,) + l-(g^^Fo,Fo,), (3.50) 

where e is such that < e < 27. On the other hand from the definition of tqo it 
follows that 

g''Fo,Foj < TOO. (3.51) 

and using p:T7|) g'^ FoiFoj < Ce^^^r^ ^hus we deduce from ^^M^ and (jXSDl) the 
following 



d 

d7 



(g'^V^Vj) < (-27 + e + Ce-^^)g''V,V, + Ce"^". 



Setting V := e^^T-e)^^,^ y.y^., it follows that 

dV 

— < Ce-^'-y + Ce-^^ 
ar 

V is thus bounded by Gronwall inequality and then 

g'^VrVj < Ce(-2T+^)^ (3.52) 

Therefore g'^^ViVj is bounded. This is enough to deduce from p.48p that for m > 0, 
we have 

so that V goes to infinity as does r. The completeness of causal trajectories is then 
proved. 

3.6 Proof of Theorem 11.61 

The proof is based on a bootstrap argument. 
By hypothesis 

KoA(to) ~ 1| < <5, |(e-V')(^o)| < <5, \{e^E){h)\ < S 
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|^^2g2;x(to) „ 3 _ 3fce2''(*o)| < s, w{to) < 5, 
and by continuity, this implies that 

\t\{t) - 1| < 26, \ie-^^i'){t)\ < 26, Ke^E){t)\ < 26 
|^^2g2M(t) _ 3 _ 3/je2'^(*)| < 26, wit) < 26, 

for t close to to. 

Let Ci and e be constants for < Ci < 1 and < e < 1/2. We can reduce 6 
if necessary so that 26 < Cit^^^'^. Then there exists some time interval on which 
the following bootstrap assumption is satisfied 

\tX - 1| < \e-^^l'\ < Cit-^+\ \e^E\ < Cit-^+' (3.53) 

lAi^e^^ - 3 - 3fce2^| < Cit-^+\ w{t) < Cit-^+\ (3.54) 

Consider the maximal interval [to,t*) on which (j3.53p - (|3.54p hold and suppose 
is finite. 

Let us continue with the following set of equations: 

e"^^(2iA + l) + fc- At^ = 87rt2p (3.55) 

e-^^i' = -A-Kte^'j (3.56) 

dtit^e^E) = -t'^e'^b (3.57) 

e-2''(2i/i- 1) - fc + Ai^ = BTTt^p (3.58) 

w = -Xw - e^^-^fi' y^l + + L/t^ T e^+'^E. (3.59) 
From (|3.55p we have 

a - 1 = i (Ae^^t^ _ 3 _ 3;,g2M) ^ ^g2M ^ 47rt2e2M^^ (3 gQ^ 

from ((338)) 

5t[-ite-2M(Ai2e2M _ 3 _ 3fcg2M)] ^ _8^t2p, (3.61) 

and from p.59[) 

dt{tw) = -t-i(U - \){tw) ~ te^'{e-^n')^/l + w2 + L/t2 zp te^ie^E). (3.62) 

Consider a solution of the full system (|1.12p - (|1.25p on the interval [to,t^) which 
satisfies the bootstrap assumption (|3.53p - (|3.54p . Putting inequalities (|3.53p - (|3.54p 
into equations p.56p - (|3.57p . (j3.60p - (|3.62p allows new estimates to be derived. For 
this purpose it is important to have estimates for the matter quantities j, p and b. 
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Let F := max{|l /+ |1 / ^^nd L the maximum value of L over the 

support of /+ or 



/OO /'OO 
w{f++ndLdw\ 



2 



< 47rLi^Ci2r'^+2% (3.63) 



< (4^iFCi2 + iCi^ + ^c)^-4+'^ (3.64) 



/OO /'OO 
/ ={f+-f-)dLdw\ 

< 47rZFCii~*+2^ (3.65) 

An estimate for e^f^ is also required. 

e^^ = A-4-2(A<2e2M) 

< A-4-2[(Ai2e2M _ 3 _ 3^62/^) + 3 + ike^''] 

< A-H-2[c^t-3+£ ^ 3 ^ 3fce2'^]. 

If A: < the last term in the latter inequality can be discarded. If fc > then we 
need 3A^^i^2 <; i [ a<2 > 3. Assume for the moment that fc < 0. Then 

< A-4-2[C^t-3+^ + 3] 

< 3A-h-^ + CiK-H-^+'. (3.66) 

It follows that 

|47rte''j| < 47ri[3A-4-2 + CiK-^r'^+'f'^ x ATTLFClt-'^+^' 

< 167r2LFC2t-3+2.[V^i-l +cl/2^-l/2^-5/2+e/2i 

keeping the worst powers and using ()3.56|) implies 



|e /I 



^■■'1 < 167r2LFC2(\/3]FT+ Cy'A-i/2)t-4+2^ C2r4+2^ (3.67) 
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Note that here there is no dependence on the initial data except for the conserved 
quantities L and F. Moreover C2 = 0{Cl). 
From (|3.65p - (|3.66p we have 

\dt{t\^E)\ = \t^e^^h\ 
keeping the worst powers gives 



and integrating this in time yields 



4+2f 



< 



=: Cst-^. 



(3.68) 



The constant C3 in the last inequality depends in a transparent way on the initial 
data. We have 

STTt^p < AiriSirLFCf + Cf + c)t"2+2% 
so that using p.6ip and integration gives 

I - ite-2^(Ai2e2^ - 3 - 3fce2'^)| < itoe"^^^*"^ lAt^e^'^^*") - 3 - Sfce^'^^*")] 
3 3 

+ l^^^i^-LFCf + Cf + c)t-'+'^\ 
At this point the assumption e < 1/2 is needed. Then 

|-iie~2A.(^^2g2M_3_3fcg2M)| < h„\Atl-3e'^'''-''>^-3k\ + -^{8TTLFCf+Cf+c). 
Using p.66p and keeping the worst powers, it follows that 



< [3A-1 + CiA-i] 



I2n 



to\Atl - 3e"2M(*«) - 3fc| + ——{8T:LFCf + Cl + c) 



(3.69) 
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Now let us examine the evolution of w. 

Using (|3.66p . (|3.74[) . (|3.68[) . the bootstrap assumption, the fact that Ci < 1 and 
then keeping the worst powers gives 



< 



and using (|3.62p . the bootstrap assumption and integration it follows that 
tw{t) < tow{to) + I Cis~^+^sw{s)ds 

J to 



J to I 1 - e 

< tow{to) + 6tge-^(*°)|£;(io)| + SdnLFA-^^^Cf + 87rA-i/2(2 + VL)^ 

+ ^ =: C,. (3.70) 
An estimate on the matter quantity p is needed. Using (I3.70p gives 



00 />00 



y/l + w^+L/t^{f+ + f-)dLdw 



00 Jo 



< 

< 

< 



t'^ J -w JQ 

A-kLF 



L/mFdLdw 



AttLF 



w 



C'st-^Jl + Clt-^ + Lt-^ 



< A-KLFC:-,t-'\l + C^t 



(3.71) 



< AnLFCrJ-^il + Cst'^ + 



Lt-') + \clt-^+'^ + \ct- 



(3.72) 



31 



In the case of plane symmetry fc = equations (|3.66p . p.69p and p.72p imply that 

1 /I -LT^ 1 



< 



+ mTT^LFA-^diS + C7i)(l + ^5 + + ^Cf + ic 



(3.73) 



Whereas in the case of hyperbolic symmetry k = —1, A ^{3 + Ci)t ^ appears as 
an adding term in the right hand side of equation (I3.73p . i.e. 



< 



^Ci + 16n^LFA-^C5i3 + Ci)(l + Cg + /E) + ^C^ + ^c + A'^S + Ci) 

(3.74) 



For the latter inequalities we only kept the worst powers. We have 

|^A-1|<C6^~^ (3.75) 

where Cg is the constant written out in p.73p for the case fc = 0, and in (|3.74p for 
the case fc = — 1. 

The constants C2-Cq appearing in equations (|3.70p - (|3.75p are all less than or 
equal to C x {g{6) + Cf), with C a positive constant and g{S) a positive function 
of 6 tending to as 5 tends to 0. Therefore it is always possible to choose Ci 
and 6 small enough in such a way that CCi < 1/2 and Cg{6) < Ci/2, and so 
the constants C2-C% are all less than Ci. This closes the bootstrap argument as it 
implies that (|3.53p - (|3.54p hold on an interval [to, ^i), with ti > t*. This contradicts 
the maximality of the interval [t{j,t<t). Therefore — oo. To complete the proof 
of theorem 11.61 it remains to show that the spacetime is complete. In fact recall 
that as in |12| the relation between coordinate time t and proper time t along the 
trajectory is given by 

dr 



The decay estimate on e^E and equations (|2.8p and (|1.2ip can be used to obtain 
the inequality 

e^'" > ^ T— , for fc < 0. 

~ C + {C-k)t + f i3 

It follows that 

> cr^, t > to- 
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Thus 

dr ^ Ct'^ 
'dt - y/m^ + C + L' 
and so r goes to infinity as does t. Theorem II .61 is then proved. 
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